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Abstract. We consider the SO (3) lattice gauge theory at weak coupling, in 
the Villain action. We exhibit an analytic path in coupling space showing 
the equivalence of the SO(3) theory with SU{2) summed over all twist 
sectors. This clarifies the "mysterious phase" of 50(3). As order parameter, 
we consider the dual string tension or center vortex free energy, which we 
measure in SO (3) using multicanonical Monte Carlo. This allows us to 
set the scale, indicating that 0(7OO) 4 lattices are necessary to probe the 
confined phase. We consider the relevance of our findings for confinement 
in other gauge groups with trivial center. 



1. Motivation 

Our main motivation for a numerical study of the 50(3) gauge theory on 
the lattice is to clarify an apparent paradox. 

On one hand, 50(3) = SU(2)/Z2, so that the 50(3) theory is an 
SU(2) theory in the adjoint representation. On the lattice, the Yang-Mills 
action —\f d A xF^ can be obtained as the naive continuum limit of a 
plaquette action taken in any representation. In the usual Wilson action 
PfJ2p(1 ~ ^^fUp), the trace of the plaquette matrix Up is taken in 
the fundamental representation. The 50(3) theory corresponds to action 
Pa2~2p(1 ~ \ r ^AUp), with the trace taken in the adjoint representation. 
The universality of the continuum fixed point leads us to believe that the 
50(3) and SU(2) theories are equivalent not just in the naive continuum 
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limit, but also non-perturbatively. Thus, 50(3) should confine (at low tem- 
perature T < T c ), just like 5*7(2). 

On the other hand, Center (SO (3)) = 1. There is no center symmetry to 
break in the case of 50(3). This means that the well-accepted connection 
between center symmetry breaking and the deconfinement transition [1] 
cannot apply. Another order parameter must be found. In the process, we 
may learn some general lesson about the gauge group structure necessary 
to sustain confinement. 

Our attention will focus on topological excitations common to both 
5*7(2) and 50(3): center vortices. Remarkably, a group with a trivial cen- 
ter may still support center vortices. A center vortex is a two-dimensional 
topological excitation. To prevent its action from diverging, it is a pure 
gauge at r = oo, characterized by the gauge transformation g{0) applied 
to the trivial vacuum. In SU(N), the non-trivial topology comes from the 
possibility that g(0 + 2ir) = g(9)exp(i^k),k = 0,..,N - 1. The integer 
k is the "twist", defined mod TV, and the Wilson loop at r = oo takes 
value W(r = oo) = g- 1 (0)g(2ir) = exp(i^fe). Notice that both W(r = oo) 
and the gauge field A^r = oo,6) = g [9)d^g(9) are unchanged if g(9) 
is multiplied by a Zn center element. This implies that our topological 
excitations are characterized by equivalence classes SU(N)/Z^ of gauge 
transformations. Twist arises from the non-trivial mapping of the r = oo 
circle to these equivalence classes: Tli(SU (N) / Z n) = Zjy. For a general 
non-Abelian gauge group G, center vortices exist if Hi(G/C enter (G)) ^ 1. 
Thus, both 5*7(2) and 50(3) admit center vortices. 

One may also realize this fact by recalling that center vortices arise as 
local gauge singularities after center gauge fixing [2], where center gauge is 
just the Landau gauge in the adjoint representation: the gauge condition 
only makes use of the 50(3) content of the gauge links, so that the gauge 
singularities only depend on the 50(3) part. 

2. Center vortices in 5*7(2) 

The relevance of center vortices for confinement in 5*7 (2) has been studied 
numerically in two ways. One can fix the gauge to a maximal or Laplacian 
center gauge, where center vortices can be identified as gauge singularities 
[2] or as P-vortices [4] and center projection can be performed. One can 
also study center vortices in a gauge- invariant way via the 't Hooft loop [3] , 
since a 't Hooft loop of contour C creates a fluctuating center vortex sheet 
bounded by C [5] . In the latter 't Hooft loop of maximal size L p x Lo- 

is equivalent to enforcing twisted boundary conditions [3] in the orthogonal 
[iv plane. A simple way to see this is to observe that the twist matrices 
IL^rjj, satisfy (for 5*7(2)) fi^fi,, = — f^lL^. Therefore, a Wilson loop of 
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Figure 1. The effect of temporal twist in SU(2): ratio Z t bc/Z p bc as a function of tem- 
perature, for various lattice sizes (from [6]). 

size x L v picks up a factor Q, p £l u Q)^l\, = —1. In other words, twisted 
b.c. multiply the largest Wilson loop at "r = oo" by a center element, 
precisely as a center vortex does: twisted b.c. (tbc) create one center vortex 
(mod 2) as compared to periodic b.c. (pbc). 

The ratio of partition functions Z t b c /Z p bc = exp(— Fqv) probes the free 
energy Fqv of a center vortex. Note that UV divergences cancel out in the 
ratio, so that it is a well-defined continuum quantity. On a hypercubic L 
lattice, this ratio rapidly tends to 1 as L increases above ~ 0.8 fm [7]. At 
finite temperature, considering electric twist (i.e. in a temporal plane), this 
ratio is an order parameter for confinement [6]. Below T c , the correlation 
length £ is finite and boundary conditions do not matter if the system size 
L is S> £: Ztbc/Zpbc — > 1. Above T c , center vortices are "squeezed" by the 
compact Euclidean time. Their free energy increases as a(T)L 2 , generating 
an area law for the spatial 't Hooft loop. The prefactor cr(T) is the dual 
string tension, which varies with T and vanishes as T — > T c + with Ising-like 
critical exponents [8]. Since we claim to measure in Ztb c /Z p b c a physical 
quantity, it should be measurable as well in the SO (3) theory. We will see 
that this is indeed the case. Z t b c /Z p b c is an order parameter also for SO (3), 
and cr(T) provides a way to measure the temperature also in that theory. 

3. SU{2) — SO(3) phase diagram and the 50(3) "mysterious phase" 

Contrary to the Wilson action /3p X)p(l — \ ^ t fUp) which enjoys a smooth 
crossover from strong to weak coupling, the adjoint action /?A2~Zp(l — 
^TtaUp) produces a bulk, first-order phase transition at Pa ~ 2.5, which 
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prevents an arbitrary decrease of (3a while studying the weak coupling 
phase relevant for the continuum limit. The common wisdom is therefore 
that 50(3) gives the same continuum physics as SU(2), but that the lat- 
tice spacing is kept small by the necessity of staying on the weak coupling 
side of the bulk transition. One purpose of our study is to check this sce- 
nario. There is room for skepticism, because the phase diagram of a mixed 
(Pf,Pa) action, first studied by Bhanot and Creutz [9], shows the weak 
coupling (3p = (i.e. 50(3)) phase qualitatively separated from the SU(2) 
theory by a line of first-order phase transitions: there is no analytic (/3p, (3a) 
path connecting the two. The possibility that the two theories are different 
has actually been proposed, with numerical results to support it, in [10]. 

Another feature of SO(3) encouraging such speculation is the existence 
of a "mysterious phase", reported by Datta and Gavai [11]. This phase ap- 
pears in all respects similar to the "normal" phase, except for the Polyakov 
loop P: instead of ^TrpP approaching ±1 as (3a — » oo, it approaches zero 
(so that, in the adjoint representation, Tr^P —> — 1 instead of +3). Tunnel- 
ing between the two phases appears extremely infrequent, increasingly so 
at higher (3a- One wonders how the SO (3) theory could have two distinct 
continuum limits, both similar to SU(2). 

To investigate this puzzle, we consider the Villain action for the SO (3) 
theory, which associates a discrete variable ap = ±1 to each plaquette 
and has action S = —(3vJ2p a P^ T FUp. Note that ap can be integrated 
analytically, giving S = — J2p l°g cosh(^TrpUp), which is manifestly in- 
variant under a change of sign of Tr pUp, like the usual adjoint action 
PaHp^-Y^aUp) = |/?AEp(l-(iTr F f/ P ) 2 ). Therefore, we deal with a 
genuine 50(3) action, which shows the same bulk first-order transition (at 
(3y ~ 4.47) and the same "mysterious phase" [11]. The advantage of the 
Villain choice is that the connection 50(3) <-> SU(2) is easier to display. 

4. Solving the puzzle: an analytic path between SU{2) and 50(3) 

This section contains no new results. We are simply pulling together obser- 
vations previously made by various authors. Even the analytic path between 
50(3) and 51/(2) was already hinted at by Halliday and Schwimmer [12]. 

The first observation, originally due to Mack and Petkova [13], was 
made precise by Tomboulis and Kovacs [14]. The SU(2) partition function 
for the Wilson action can be exactly rewritten as an 50(3) partition func- 
tion in Villain form, with the additional constraint that ap-monopoles are 
forbidden. Namely: 



Zsum = j VUePXpW' = c £ J VU^p a ^ u "HS(llap-l) 

ap=±l cubes 6 

(1) 
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where c comes from the normalization of the ap integration. 

The second observation is due to Alexandra and Haymaker [15], and is 
a refinement upon the first. While eq.(l) applies in an infinite volume, on 
a 4-torus one needs some additional global constraints on the ap variables: 

N^= J] «P = +1 (2) 

P 6 plane \iv 

for each product of x L u variables ap in each \iv plane. Only configura- 
tions satisfying (2) can be mapped to SU(2) configurations with periodic 
b.c. Note that, in the absence of ap-monopoles, ]lp e plane ^v a P ls the 
same for every parallel plane [iv. Therefore, eq.(2) defines 6 constraints 
only. If one imposes = — 1 for some [iv orientation, then the mapping 
(1) is preserved, except that the SU(2) b.c. are twisted in the [iv plane. 

It is then straightforward to remove the constraints (2) by summing 
over all 2 6 N^ v possibilities. One obtains: 

£ / VUe^ P ^ F u P -Q 6([[ap _ 1) = c -i £ Zsum (3) 

ap=±l cubes 6 all twist sectors 

The SO (3) Villain partition function, with monopoles removed, is just the 
sum of 577(2) partition functions over all possible twisted b.c. 

The "mysterious phase" then must be one (or more) of these twisted 
SU (2) sectors. Indeed, one knows how to construct a so-called "twist-eater" , 
a groundstate with zero action and twisted b.c. [16]. Suppose N xt = —1 is 
enforced through ap = —1 for all xt plaquettes at (1, y, z, 1), +1 otherwise. 
Then, set all links to 1, except U x (l,y, z, t) = io\ and Ut(x,y, z,l) = iui- 
This choice yields ^Tr pUp = +1, except for x = t = 1 in xt planes, where 
Up = —1. This minus sign cancels with that of ap to give zero action. Note 
that the two Polyakov loops P x and P± in the x and t directions are now 
io\ and zo^j satisfying TrpP = 0. 

Therefore, we guess that, whenever the 50(3) system is in the myste- 
rious phase characterized by Tr pP t ~ (i.e. TY^i^ ~ — 1), twist in the ap 
variables must be present in some of the temporal planes. This proves to 
be correct, as illustrated in Fig. 2, where we show the Monte Carlo history 
of TtAPt, accompanied by that of k Xjy>z , where 

1 / 1 LyLz \ 

is the average, mapped to the interval [0, 1], of the op-twist N xt in each 
xt plane. This averaging is necessary because our SO (3) action does not 
forbid ap-monopoles. They have a small (oc e~ mlS ) density, which allows 
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Figure 2. Monte Carlo history of the adjoint Polyakov loop (bottom) and of the 3 electric 
twist variables (top). The trace of the adjoint Polyakov loop is negative whenever twist 
is present (4 4 lattice, (3 = 4.5). 

N x t = Yip ap to fluctuate from one parallel plane to another, so that k x y z 
fluctuate near or near 1. Note also the very slow Monte Carlo dynamics, 
which prevent an ergodic sampling of configuration space on all but the 
smallest lattice sizes unless one resorts to special strategies, see sec. 5. 

By measuring Polyakov loops in other directions, one can resolve whether 
one or more of the k X y z are near 1, and obtain a one-to-one mapping be- 
tween Polyakov loops and twist sectors [17]. 

Completing an analytic path between 50(3) and SU(2) requires one 
last, simple step. We need to allow again ap-monopoles, which can be done 
via a chemical potential A, giving the partition function 



In the limit A — > +oo, ap-monopoles are forbidden and eq.(3) holds. The 
SO(3) Villain action which we study corresponds to A = 0. The phase 
diagram in the (A,/3) coupling plane has been studied in [11]. We repro- 
duce it in Fig. 3. At strong coupling, monopoles condense: their density 
is 0(1), which introduces strong fluctuations in Jlp^-P from one parallel 
plane to another. Variables k x ^ y>z fluctuate around |, and the concept of 
twist is ill-defined. On the other hand, at weak coupling, monopoles are ex- 
ponentially suppressed, long-range order appears among the ap variables, 
and twist sectors become well-defined. There is no non-analyticity between 
(A = 0, (3 > (3 C « 4.47), which is the weak-coupling SO (3) Villain theory un- 
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Figure 3. Phase diagram in the coupling plane (A, /3), where A is the monopole chemical 
potential and /3 the 5*0(3) gauge coupling (adapted from Ref.[ll]). 



der study, and A = +00, which is equivalent to SU (2) summed over all twist 
sectors as per eq.(3). Therefore, the non-perturbative continuum physics of 
SO(3) is the same as that of SU(2), summed over all twist sectors. 

What remains to do is to put this statement on a quantitative footing, 
by matching the lattice spacings o(/3) in 50(3) and SU(2). We accomplish 
this next by measuring the dual string tension <r(T) in both theories. 

5. The order parameter and its measurement 

Defining an order parameter for confinement in SO (3) is not a trivial mat- 
ter, as emphasized by Smilga [18]. The trace of the SU(2) Polyakov loop, 
like that of all Wilson loops, is identically zero in the fundamental repre- 
sentation. In the adjoint representation, the trace of the Polyakov loop is 
always non-zero, since adjoint charges are screened by gluons at any tem- 
perature. Here, we use the dual string tension o"(T), which is zero in the 
confining phase T < T c , positive and increasing with T in the deconfined 
phase T > T c , and which we already showed to be a good order parameter 
in the SU(2) theory (see [8, 6] and Fig. 1). 

In SU(2), on a lattice of physical size L 3 x T , the dual string tension 
can be obtained from the ratio of twisted (in a temporal L x T _1 plane) 
and periodic partition functions: 

a(T) = lim -(1/L 2 ) log Z tbc /Z pbc (6) 

L^oo 

The presence or absence of twist is imposed by hand, by a specific choice 
of boundary conditions. In 50(3), twist sectors are automatically summed 



Figure 4- The three-dimensional reweighting table used in our multicanonical Monte 
Carlo. No twist is at bottom left, twist in all 3 planes at top right. The table enhances the 
probability of sampling the saddles between twist sectors. Its entries vary by 13 orders 
of magnitude (8 3 x 4 lattice, /3 = 4.5). 

over, and twist becomes an observable, similar to the usual topological 
charge in SU(2). Sectors of different twist, which would be disjoint in the 
continuum limit, are connected at finite f3 via saddle points which are lat- 
tice artifacts, analogous to topological "dislocations". The assignment of an 
integer twist (in each plane) to a given configuration is subject to ambigui- 
ties similar to the assignment of a topological charge. Fortunately, we work 
at a very weak coupling, so that such ambiguities have a negligible impact 
in our case: the twist variables k xyz (eq.(4)) stay near or 1. The price 
we have to pay is that the twist sectors are separated by very high action 
barriers AS, nearly impassable for an ordinary Monte Carlo algorithm. 

Our strategy to solve this technical problem is multicanonical sampling 
[19]. Instead of assigning probability \eT s to each configuration, we multi- 
ply this probability by a factor T(k) which favors the saddle points between 
twist sectors k ~ and k ~ 1, tuned such that the resulting sampling prob- 
ability is flat from one sector to the next. The Monte Carlo process then 
amounts to a free diffusion across twist sectors, with dynamics accelerated 
exponentially (by a factor ~ e@^ s ). 

In our case, we want to facilitate the Monte Carlo sampling of twist in 
all 3 temporal planes. This leads us to a multicanonical sampling with a 3- 
dimensional reweighting factor T(k x , k y , k z ) which can be represented by a 
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density of states 
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Figure 5. Cut of the three-dimensional reweighting table along its diagonal (blue *). The 
result of multicanonical Monte Carlo sampling is a nearly flat histogram (red +). The 
density of states shows the dominant twist-0 and the smaller twist-3 sectors (green x). 



table since k x , y , z = k/L 2 , k = [0, .., L 2 ]. This table is constructed iteratively, 
using converged values on small lattices to form starting values on larger 
ones. Such a table is displayed in Fig. 4, for an 8 3 x 4 lattice at j3 = 4.5. 
Although (5 has been chosen as small as possible (the bulk transition to the 
strongly coupled phase occurs at [5 C sa 4.47), the necessary enhancement of 
the saddle points reaches 10 13 , as seen more clearly in Fig. 5, which is a ld- 
cut of the same table along its diagonal (from k x>y>z = (0,0,0) to (1, 1, 1)). 
There, the resulting flatness of the Monte Carlo sampling probability is 
clearly visible. The measured density of states shows two well-separated 
peaks, corresponding to twist sectors (0,0,0) (analogous to SU(2) pbc) and 
(1, 1, 1) (analogous to SU(2) the in all 3 temporal planes). The strong sup- 
pression of the saddle point confirms that ordinary Monte Carlo sampling 
would remain hopelessly "stuck" in one sector, as observed in earlier stud- 
ies. However, in spite of the great multicanonical acceleration, the Monte 
Carlo evolution of the twist variables k X)V ^ z is still slow, and simulating a 
10 3 x 4 lattice remains beyond the edge of our computer resources. 

This prevents us from reaching spatial sizes large enough for a reliable 
extrapolation of eq.(6) to L — > co. Our extrapolation depends on the ansatz 
we choose for the finite-size effects. Nevertheless, we can still calibrate the 
SO(3) lattice spacing a((3 SO (3) = 4.5), by comparing our L = 4, 6,8 results 
with those obtained in SU (2) on the same lattice sizes L 3 x 4. In both theo- 
ries, our observables are Z\ hc jZ v b c = exp(—F^ v ), where k = 1, 2, 3 refers to 
twisted b.c. in 1, 2 or 3 temporal planes. In SO(3), these observables are ob- 
tained from the Monte Carlo distribution of the twist variables k XtV:Z . This 
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Figure 6. Comparison of the twist free energies between SO(3) at /3 = 4.5 and SU(2) 
at various /3's. Electric twist in 1, 2 and 3 planes is considered, on lattices of size 4 3 ,6 
and 8 3 x 4. One finds f3 so(3) = 4.5 < — > /3 SU(2) = 4.10(15). 



distribution falls into 2 3 distinct, well-separated peaks, as in Fig. 5, and we 
measure the ratio of the population around each peak to that around the 
pbc peak at kx y Z ~ 0. In SU{2), these ratios are obtained from different 
Monte Carlo simulations, following the method of [8]. 

The twisted free energies Fq V are continuum quantities, which depend 
upon the spatial size L and the temperature T. If the 50(3) and SU(2) 
theories represent the same continuum physics, then to each /?so(3) should 
correspond a value Psu(2) j which yields the same lattice spacing and thereby 
the same Fq V for equal lattice sizes, modulo small lattice artifacts. We test 
this statement in Fig. 6. The 9 free energies Fjj V ,k = 1,2,3 measured on 
4 3 , 6 3 and 8 3 x 4 lattices in SO(3) at PsO(3) = 4-5 are compared with similar 
quantities measured in SU (2) at different values of flsucz)- The straight lines 
show a linear interpolation in Psu(2) °f the SU{2) data. One observes a very 
good match of all 9 observables, for (3su(2) ~ 4.10(15). 

Thus, a -1 (/? S o(3) = 4.5) w a~ 1 (/3 5[/{2) = 4.1), which is about 200 GeV! 
As conventional wisdom asserts, our SO(3) lattice is very fine indeed. To 
reach low temperatures T < T c and probe the confined phase would require 
a lattice of size O(700) 4 , far beyond what is currently achievable. 

Note the closeness of the matched SO(3)-SU(2) bare couplings (-y = 
4.5 vs 4.1). This should come as no surprise. Lattice perturbation theory is 
identical between the SU{2) Wilson action and the SO(3) Villain action: 
the difference resides in the ap-monopoles, which do not appear in the 
perturbative expansion. Therefore, Ai a ttice is the same in both theories, 
and one should expect similar values for the non-perturbatively matched 
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/3's. The ap-monopoles disorder the 50(3) theory slightly, which is why 
the matching SU (2) (3 is slightly smaller. 

6. Conclusion 

To summarize, conventional wisdom prevails and there is no mystery. The 
50(3) lattice theory at weak coupling gives the same non-perturbative 
physics as the SU(2) theory and a common confinement order parameter 
exists for both: the center vortex - or twist - free energy. The difference is 
that a definite twist sector is selected via the choice of boundary conditions 
in SU(2), whereas all sectors are automatically summed over in 50(3). 

Unfortunately, the numerical study of the 50(3) theory seems limited 
to very fine lattice spacings, because a bulk transition prevents exploring 
continuum physics at stronger bare couplings. One way out is to suppress 
the formation of monopoles whose condensation triggers the bulk transition 
[20]. However, this appealing approach suffers from an unpleasant, techni- 
cal side-effect. The Monte Carlo sampling over the various twist sectors, 
which is essential to measure the confinement order parameter, is made 
more difficult by higher action barriers separating the sectors. This situa- 
tion is analogous to simulations of overlap or domain-wall fermions. There, 
dislocations separating topological sectors cause trouble in the inversion of 
the Dirac operator. They can be suppressed by choosing a gauge action 
which assigns them a large action. The price to pay is a much slower Monte 
Carlo evolution of the topological charge. One may hope that a clever en- 
hancement of the sampling probability in the neighborhood of the saddle 
point, and only there, as in [21], can solve both types of problems. 

Finally, one can try to generalize the lesson learnt here about the struc- 
ture of the gauge group G necessary for confinement. 50(3) shows that 
a non-trivial center is not required. On the other hand, the existence of 
the dual string tension arises from that of the twist sectors. Those in turn 
follow from the non-trivial first homotopy group H\(G /C enter (G)), which 
is the same for SU{2) and 50(3). A non-trivial Ui(G /C enter (G)), or more 
precisely Y[\{G/Zq), where Zq is the discrete part of the center of G, is 
necessary for a dual string tension to be defined. The non-trivial elements 
of this group are the center vortex (or twist) excitations. Therefore, we 
conjecture that the existence of center vortices is a necessary condition for 
the existence of an ordinary string tension, i.e. for confinement. It is not a 
sufficient condition: compact U(l) has twist sectors but its weak-coupling 
phase does not confine. On the other hand, there are non-Abelian Lie groups 
which do not meet this condition, namely G2, F4 and Eg. Therefore, we 
conjecture that these cannot confine (in the sense that the Wilson loop 
cannot obey an area law). 
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